Critical behavior of 0(n)— symmetric systems with reversible mode— coupling terms: 

Stability against detailed— balance violation 
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We investigate nonequilibrium critical properties of C'(n)-symmetric models with reversible 
mode-coupling terms. Specifically, a variant of the model of Sasvari, Schwabl, and Szepfalusy 
is studied, where violation of detailed balance is incorporated by allowing the order parameter and 
the dynamically coupled conserved quantities to be governed by heat baths of different tempera- 
tures Ts and Tm, respectively. Dynamic perturbation theory and the field-theoretic renormalization 
group are applied to one-loop order, and yield two new fixed points in addition to the equilibrium 
ones. The first one corresponds to = Ts /Tm = oo and leads to model A critical behavior for 
the order parameter and to anomalous noise correlations for the generalized angular momenta; the 
second one is at = and is characterized by mean-field behavior of the conserved quantities, by 
a dynamic exponent z = d/2 equal to that of the equilibrium SSS model, and by modified static 
critical exponents. However, both these new fixed points are unstable, and upon approaching the 
critical point detailed balance is restored, and the equilibrium static and dynamic critical properties 
are recovered. 

PACS numbers; 05.70.Ln, 64.60. Ak, 64.60.Ht 



I. INTRODUCTION 

Nonequilibrium steady states (NESS) have been much 
investigated, the main goal being the discovery of their 
common and distinguishing features as compared to equi- 
Ubrium states. A promising approach to this problem is 
the study of phase transitions: Since equihbrium criti- 
cal phenomena display a large degree of universaUty, it is 
natural to ask to what extent these universal features re- 
main characteristic of nonequilibrium phase transitions. 

The basic complication with NESS is that, in addition 
to the interactions which entirely define the equilibrium 
properties, the dynamics is also essential in determining 
the steady state properties. Thus, for example, a classi- 
fication of nonequilibrium phase transitions requires not 
only the understanding of the role of symmetries of the 
order parameter, the range of interactions and the di- 
mensionality of the system, but the clarification of both 
the relevance of conservation laws imposed by dynami- 
cal symmetries and the range of the dynamical processes. 
Possibly, new dimensionality and anisotropy effects in the 
dynamics may also be important. 

The most frequently studied models with nonequilib- 
rium phase transitions are generalizations of systems 
with model A type dynamics The transitions in 

these systems have been shown to be robust against lo- 
cal nonequilibrium perturbations which do not conserve 
the order parameter |^ and, remarkably, this robustness 
was found to persist even if the dynamical perturbations 
broke the discrete symmetry of the system ||]. Both lo- 
cality and the nonconserving character of the perturba- 



tions are essential for the phase transition to stay in the 
Ising universality class. Indeed, nonlocal nonequilibrium 
dynamics generates effective long-range forces and thus 
changes the universality class dramatically [Q. 

The nonequilibrium generalizations of model B type 
dynamics (with conserved order parameter) are more in- 
teresting. External fields or local, anisotropic, nonequi- 
librium perturbations may drive the system into a NESS 
with phase transitions which are not characterized by 
any known equilibrium universality class or belong 
to universality classes with long-range interactions [^,0 . 

Nonequilibrium generalizations of the case when a non- 
conserved order parameter is coupled to a conserved 
quantity have been considered in Ref. H] where it was 
found that linear coupling to a conserved quantity gen- 
erates power-law correlations for the order parameter. 
This suggests that, in this situation, long-range effective 
interactions are generated in the system, which in turn 
govern the critical behavior at the phase transition. 

There are several other nonequilibrium phase transi- 
tions which have been studied without considering any 
equilibrium context. Most notable among these are phase 
transitions associated with the presence of an absorbing 
state (directed percolation) , and the roughening tran- 
sition in surface growth and equivalent models such as 
the Kardar-Parisi-Zhang equation . 

In this paper we continue the investigation of nonequi- 
librium generalizations of models originally proposed to 
describe equilibrium critical dynamics. Our aim is to 
study an example where there is a reversible mode- 
coupling between the order parameter and another (con- 
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served) field, using the field-theoretic dynamic renormal- 
ization group (RG) A simple example of this 

type of systems is the Heisenberg model for isotropic fer- 
romagnets where precession terms introduce a coupling 
among the different spin components (model J according 
to the classification in Ref. for early RG studies of 



this model see Ref. |13 ; a comprehensive review of the 
critical dynamics of ferromagnets is given in Ref. [Q). 
However, similar to the purely relaxational dynamics of 
models A and B, the effect of a (spatially isotropic) vi- 
olation of detaile d ba lance can be removed by a simple 
rescaling (see Sec.IIB). We shall thus mainly consider a 
more complicated model which was originally introduced 
by Sasvari, Schwabl, and Szepfalusy in the context of 
structural phase transitions This SSS model con- 

sists of a non-conserved n-component order parameter 
purely dynamically coupled to the n(n — l)/2 conserved 
generalized angular momenta related to the underlying 
0{n) symmetry of the system. The n = 2 realization de- 
scribes the critical dynamics of planar ferromagnets and 
superfluid Helium 4 [^6| (for reviews regarding dynamic 
critical phenomena in superfluid Helium, see Ref. [p^), 
while the case n = 3 corresponds to the dynamics of 
isotropic antiferromagnets flq]. The SSS model, with 
its dynamic exponent z — dj2 (below the upper critical 
dimension dc = 4) thus encompasses models E and G 
(according to Ref. [Q) as special cases. 

We shall generalize the previous field-theoretic RG 
studies of the SSS model |l^j2^ to a nonequilibrium situ- 
ation by assuming that the order parameter components 
and conserved angular momenta M"^ are attached 
to heat baths of different temperatures Ts and Tm, re- 
spectively. Thus the detailed-balance condition required 
for near-equilibrium dynamics is violated and the flow 
of energy between the two heat baths ensures that the 
steady state is out of equilibrium. The introduction of 
two temperatures leads to an additional variable in the 
problem, namely the temperature ratio — Ts/Tm (of 
which no analog can be constructed for model J). By 
studying the RG flow equations to one-loop order (first 
order in e = 4 — d), we find two new fixed points cor- 
responding to the cases 6=0 and 9 = oo, respec- 
tively, in addition to the equilibrium fixed points of the 
SSS model. The latter are (i) the usual Gaussian fixed 
point (describing static and dynamic mean-field behav- 
ior, 2 = 2), (ii) the model A fixed point, correspond- 
ing to a decoupling of the conserved fields from the or- 
der parameter, with the nontrivial static exponents of 
the 0(n)-symmetric 0^ model and dynamic exponent 
z = 2 + 0(6^), and the three nontrivial SSS dynamic fixed 
points consisting of the two so-called weak-scaling fixed 
points with the order parameter and conserved quanti- 
ties fluctuating on different time scales, characterized by 
the exponents (in) Z5 = 2 - 2(n - l)e/(2n - I) + ©(e^), 
ZM = 2 - e/{2n - 1) + 0{t^), and (iv) = 2 and 
zm = d—2, and finally (v) the strong-scaling fixed point 
with zs = Zm = z = d/2. The results for (iv) and (v) 
actually hold to all orders in e, and follow from the exact 



sum rule zs + zm = d [|l6|,yjj2l| (see Sec. pi A| ). Stability 
analysis shows that to one-loop order only the strong- 
scaling fixed point (v) is stable; however, at least for 
n ~ 2 the actual fixed-point values are rather close to its 
stability boundary, which allows for the possibility that in 
fact for superfluid Helium strong scaling may be violated 
at the Lambda transition |l^] (a two-loop study of model 
F, combined with Borel-resummation techniques, actu- 
^ly suggests the stability of a weak-scaling fixed point 

The stability of the above fixed points may change in 
a nonequilibrium situation where the order parameter 
and conserved variables are allowed to fluctuate at dif- 
ferent temperatures, which explicitly introduces different 
characteristic time scales. Indeed, while one of the two 
new nonequilibrium fixed points, corresponding to (a) 
Q — 00, is described by model A dynamics zs = zm = 2 
(with the usual 0^ model statics), albeit accompanied 
by anomalous noise correlations for the conserved fields^ 
the second new fixed point, characerized by (b) 8 = 0, 
yields, actually to all orders in e, zs = d/2 for the order 
parameter as in equilibrium, but zm = 2, i.e.: ordinary 
diffusion for the generalized angular momenta (note that 
the above-mentioned equilibrium sum rule does not hold 
here); this unusual behavior is supplemented by anoma- 
lous order parameter noise correlations^ and even modi- 
fied static critical exponents. However, stability analysis 
reveals that in fact both these fixed points (a) and (b) 
are unstable, and for any initial value of < < 00 
the fiow asymptotically leads to the stable strong-scaling 
equilibrium fixed point of the SSS model (see Sec. HI B| ) . 
Thus we conclude that while violation of detailed balance 
might be conceived as a relevant perturbation, in fact the 
underlying 0{n) symmetry in conjunction with spatial 
isotropy and the growing correlation length as the phase 
transition is approached, effectively restores detailed bal- 
ance (described by the fixed point with = 1), and thus 
asymptotically yield the usual static and dynamic critical 
behavior of the equilibrium SSS model. 

This paper is organized as follows. In the following 
Sec. H, we briefly review the derivation of Langevin equa- 
tions describing the critial dynamics of 0(n)-symmetric 
models including reversible mode-coupling terms, and 
consider the possible relevance of detailed-balance vio- 
lation for the relaxational models A and B, as well as 
model J and the SSS model. Sec. HI will then be de- 



voted to the RG study of the nonequilibrium SSS model 
as outlined above, starting with stating some general ex- 
act relations and Ward identities, followed by a detailed 
study of the one-loop perturbation theory, the ensuing 
flow equations, and a discussion of the physical content 
and stability of the RG fixed points. Finally, in Sec. IV 
we shall summarize our results again, draw some conclu- 
sions and provide an outlook on possible future research 
along the path followed in this paper. In the Appendix, 
we provide a list of the explicit results to one-loop order 
for the two-, three-, and four-point functions required 
for the renormalization of the nonequilibrium SSS model. 
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II. CRITICAL DYNAMICS OF MODELS WITH 
REVERSIBLE MODE-COUPLING TERMS 

A. General considerations 

The universal static critical behavior of a system which 
is invariant with respect to rotations of its n-component 
order parameter and displays a second-order phase tran- 
sition is described by the following 0(n)-symmetric 
Landau-Ginzburg- Wilson hamiltonian in d space dimen- 
sions 



H[{So}] 



{n T 
a— 1 Q— 1 

n 



4! 



(2.1) 



here ro = {T — T^)/T^ denotes the relative distance from 
the mean-field critical temperature T^, and we denote 
unrenormalized quantities by a subscript "0" . This effec- 
tive free energy determines the equilibrium probability 
distribution for the vector order parameter Sq, 



-Hl{SS}]/ksT 



(2.2) 



and furthermore provides the starting point for the con- 
struction of the field-theoretic static renormalization 
group which by virtue of a perturbation (loop) expan- 
sion in the nonlinearity uo provides a systematic means 
to compute the two independent static critical exponents 
r] and v either in an e expansion about the upper critical 
dimension dc = 4, or directly in fixed dimensionality d 
p3| . Here, rj is the anomalous dimension which describes 
the power-law decay of the order parameter correlations 
at the critical point, {S"{x)S'^ {x')) cx l/|x - x'I-^-^+t,^ 
and the exponent characterizes the divergence of the 
correlation length as Tc is approached, ^ oc |T — T^'" . 

The simplest dynamics that may be imposed on the 
order parameter fluctuations 6*9 (x,t) in order to de- 
scribe how the system relaxes to equilibrium (for which 
the mean-field stationarity condition (5_ff [{S'g }]/(5S'o = 
holds) is then given by the following Langevin-type equa- 
tions of motion 

-A,(.vr™+r(x,q, ,2,3, 



dt 



where the temporal average of the stochastic forces is as- 
sumed to vanish, (('"(x, i)) = 0. In equilibrium, further- 
more an Einstein relation connects the second moment 
of the uncorrelated (white) noise with the relaxation co- 
efficient, 



(r (x,t)C''(x', i')) = 2XokBT{iVrS{x - x')d{t - t')S' 



this ensures that the probability distribution PjjS"}] fi- 
nally approaches the equilibrium distribution ( |2.2|) in the 
limit t ^ 00, as can be readily checked with the aid of 
the associated Fokker-Planck equation. Eq. (2.3) incor- 
porates both the case of a nonconserved order parameter 
with purely relaxational dynamics (a = 0) and a con- 
served order parameter which as a consequence of the 
ensuing continuity equation relaxes diffusively (a = 2). 
In the classification scheme of Hohenberg and Halperin, 
these situations are referred to as models A and B, re- 
spectively, and the corresponding dynamic critical expo- 
nents describing the critical slowing down near the phase 
transition (characteristic timescales diverge as tc oc oc 
\T — Tc\^''^) are given in terms of the static exponent 77 
by z = 4 — ?7 (model B) and z = 2 -I- C77 (model A) |jl|]. In 
the latter case, however, c is a new universal number and 
therefore z is an independent exponent not determined 
by the static critical exponents. 

One may already anticip ate that an isotropic viola- 
tion of the Einstein relation (2.4), which is a consequence 
of an underlying detailed balance condition, by choosing 
a coefficient Aq instead of XoksT for the noise correla- 
tor, merely amounts to a change in the order parame- 
ter temperature T. Therefore, as long as one remains 
sufficiently close to the critical point, the universal criti- 
cal behavior (exponents, amplitude ratios, etc.) will not 
be affected, while only nonuniversal amplitudes become 
modified through a rescaled nonlinear coupling uq (see 
SeclHBland Ref. Q). 



However, in an 0(n)-symmetric system there are al- 
ways additional slow diffusive modes present. In our case 
these modes are associated with the conserved general- 
ized angular momenta Mq^ which generate the rotations 
in order parameter space. Generically, they coupl e to 
the order parameter fluctuations, and therefore Eq. (2.3) 
does not correctly describe their dynamics. Two cases 
can now be distinguished: (i) The vector order param- 
eter itself is identical with the generators of the group 
0(n); this yields, for n = 3, precisely the dynamics 
of isotropic Heisenberg ferromagnets ||l3[; model J ac- 
cording to Ref. 1^]. (ii) The order parameter is non- 
conserved, and the conserved angular momenta consti- 
tute new dynamical variables; this defines the 0{n)- 
symmetric model introduced by Sasvari, Schwabl, and 
Szepfalusy ||l^, and encompasses both model E for the 
dynamics of the XY model, i.e., of planar ferromagnets 
and superfluid Helium 4 (n = 2) [|l6| , and model G for 
isotropic antiferromagnets (n — 3) |ll8|| . 

Upon collecting the order parameter and angular mo- 
mentum components in a large vector = {S" , M"^), 
the general structure of the ensuing Langevin equations 
reads Of! 



9V"(x, t) 
dt 



Slp^'iK, t) 



(2.4) 



+ r(x,i) , 

(2.5) 
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where L" = X or L" = —D'V^ for all the nonconserved 
and conserved fields, respecti vely . The second term on 
the right-hand side of Eq. (2.5) describes irreve rsib le 
relaxation processes as in models A and B [Eq. (2.3)]; 
the first term, on the other hand, consists of reversible 
"mode-couplings", which are given entirely by the Pois- 
son brackets oc {^",'0'^}. As can be shown 

with the Kawasaki-Mori-Zwanzig projector formalism, 
assumes the form of a "streaming velocity" in 
the space of the ip" , namely 



v"[{r}] 



SQ 



af3 



(2.6) 



Note that the mode-coupling constants g are indepen- 
dent of a, which guarantees that l/"[{-(/,"}]e--f^[W°>l/'=BT 

is divergence-free. 



- (y^KV-^lJe"^'^^"^!/^^^) ^ , (2.7) 



and therefore the equilibrium distribution PcqilV""}] oc 
g-^f[{V'°}]/fcBr ig j^ot affected by the mode-coupling terms 
which are of purely dynamical origin. 

We defer the explicit construction of the mode- 
coupling terms for model J and the SSS model to the 
following subsections, and close this general discussion 
with a brief outline how one may construct an effective 
field theory from Langevin equations of the type 



9i/'"(x, t) 
dt 



-if"[{V"}](x,i)+r(x,t) , (2.8) 
and the general noise 



see Eq. (|2J), with (C"(x,i)) 
correlator 

(r(x,i)C'^(x',t')) =2L"<5(x-x')(5(t-t')^"^ , (2.9) 

see Refs. This form of the white noise may be 

inferred from a Gaussian distribution for the stochastic 
forces 



W[{C}] oc exp 



(2.10) 



eliminating via Eq. (|2.8D then immediately yields the 
desired probability distribution for the fields '0", 

w[{c}]v[{c}] = p[{r}M{r}] oc e«[w°"i?[{v^"}], 

(2.11) 

with the Onsager-Machlup functional 



G[{r}] 



x(i") 



dt 



^\ dt 

a. ^ 



K'^Hr}]) . (2.12) 



From this functional one could already construct a per- 
turbation expansion for correlation functions of the fields 
^p"; however, as for conserved quantities the inverse of 
the Onsager coefficient L" is singular, and furthermore 
high nonlinearities oc appear, it is convenient 

to introduce Martin-Siggia-Rose auxiliary fields via a 
Gaussian transformation to partially linearize the above 
functional. This leads to 



(2.13) 



with the Janssen-De Dominicis functional 



j[{r},{r}]= jd^x fdtj2 

•' a 

fdip' 



-r 



\ dt 



(2.14) 



Eq. (2.14) will provide the starting point for our dis- 



cussion of the nonequilibrium dynamics of the isotropic 
ferromagnet (model J) as well as that of the SSS model 
in the subsequent subsections. In Sec. Ill, we shall use 



the corresponding Janssen-De Dominicis functional for 
the construction of the dynamical field theory of the 
SSS model with broken detailed balance, and there- 
from infer its RG flow equatio ns. We finally remark 



that both in Eqs. ( 2.12 ) and ( 2.14 ) we have omitted 
contributions stemming from the functional determinant 
I?[{C"}]/X>[{V'"}]. As is shown in Refs. (l|jl|, these 
terms precisely cancel any acausal Feynman diagrams 
for the dynamic response function that could be con- 
structed from the above functionals; and upon restricting 
the perturbation expansion to those contributions which 
are consistent with causality requirements, we may there- 
fore safely neglect these additional terms. 



B. Model J — isotropic ferromagnets 



We now turn explicitly to the construction of the 
Langevin equation for the critical dynamics of isotropic 
ferromagnets [Q. In this case, n = 3, and the order 
parameter consists of the three spin components , S^, 
and S^. The total magnetization is a conserved quan- 
tity (hence a — 2), and in fact the S" are identical with 
the generators of the rotation group 0{n): M^^ — S^, 
Af23 ^ S"", and M^^ = -S"^. The Poisson brackets be- 
tween the spin components read 



(2.15) 



which immediately yields the streaming velocity 
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SH[{S-}] 
SSf' 



(2.16) 



for the contractions of the fuUy antisym metri c tensor 
e°'^'< with aU the symmetric terms in Eq. (|2.l|) vanish, 
leaving only the contribution stemming from the gradi- 
ent te rm in the hamiltonian. The mode-coupling terms 
(2.16) represent the spin precession in the effective field 
generated by the other spins, and in the ordered phase 
leads to propagating spin waves (Goldstone modes) with 
quadratic dispersion oj{q) oc q^. 

The complete Langevin equation for the conserved or- 
der parameter of isotropic ferromagnets (model J accord- 
ing to Ref. 01) finally reads 



dsi 

dt 



-go 



(2.17) 



with (C"(x,<)) = and 

(r(x,t)C''(x',i')) = -2Ao V25(x - x')<5(i - t')<5"'' ■ 

(2.18) 

Here we have already allowed for a violation of 
the detailed-balance condition by introducing a noise 
strength Ao that is not necessarily equal to Xok^T, where 
Ao is the spin diffusion constant. However, the form of 
Eqs. (2.17) and (2. IS) already suggests that similar to 
the case of the purely relaxational models A and B, the 
ratio Ao/Ao may be absorbed into a rescaled temperature 
T, and modified nonlinear couplings uq and go. 

This can be readily seen by employi ng th e corre- 
sponding Janssen-De Dominicis functional ( 2.14| ); for our 
nonequilibrium model J this becomes a sum of the dy- 
namic functionals for the relaxational models 



dt 



^0 



(2.19) 
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with a = 2, and the additional contribution stemming 
from the reversible spin precession term, 



Jn.c[{SS}, {S^}] = -go J d'x JdtJ2 e'^'^^S^S^V'S^o . 

Q:,/3,7 

(2.20) 

Rescaling the fields according to 



1/2 



QO. 




(2.21) 



then renders the noise strength and the relaxation con- 
stant in the quadratic part (first and second line) of 
Eq. (2.19) equal, and if in addition the rescaled static 
and dynamic nonlinear couplings 



^ Ao 
uo = — uo 



1/2 



50 



(2.22) 



are introduced, the ensuing Janssen-De Dominicis func- 
tionals for the above nonequilibrium generalizations of 
the relaxational models as well as model J appear in pre- 
cisely the same form as in equilibrium where detailed bal- 
ance holds. As both the renormalized counterparts of Uq 
and (Jo/Xq approach universal fixed-p oint v alues near the 
transition, the modifications in Eq. ( 2.22| ) merely enter 
nonuniversal amplitudes. It is therefore established that 
the critical properties of neither the relaxational models 
A and B nor isotropic ferromagnets (model J) are affected 
by violating the detailed-balance condition. It is, how- 
ever, important to note that both the 0(n) symmetry and 
the spatial isotropy of the models have been left intact by 
the above nonequilibrium generalization. For the dynam- 
ics of Heisenberg ferromagnets, we finally remark that 
the dynamic critical exponent becomes z = {d + 2 — r])/2, 
as a consequence of a Ward identity ste mming from the 
underlying 0(3) symmetry (see also Sec. [HA). Further 
details regarding the dynamic critical behavior of ferro- 
magnets may be found in Ref. [Q. 



C. The SSS model — planar ferromagnets, isotropic 
antiferromagnets 



More interesting for the issue of violating detailed bal- 
ance will clearly be a situation where there are two inde- 
pendent temperature scales conceivable, and therefore a 
simple temperature rescaling will not suffice to render the 
field theory identical to the equilibrium one. We therefore 
consider a nonequilibrium version of the 0(n)-symmetric 
SSS model, where a nonconserved n-component order 
parameter couples to n{n — l)/2 conserved generalized 
angular momenta jl^ ; possible realizations of this arc (i) 
for n — 2: the critical dynamics of the XY model llfl] 
(model E according to Ref. Q), with the order parame- 
ter components and S^, and the conserved quantity 
M^'^ = S^, which generates rotations in the xy-plane; (ii) 
for n = 3: the dynamic critical behavior of isotropic anti- 
ferromagnets, with 5*^, S^, and representing the com- 
ponents of the staggered magnetization, and M^^ = M^, 
Af23 _ j^x^ g^j^j j^,^i3 _ denoting the components 

of the magnetization itself, which are conserved and can 
be identified with the generators of 0(3) (model G |p8|). 
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The variables Mq^ are noncritical quantities, and their 
couphng to the order parameter fluctuations Sq is of 
purely dynamical character. Hence it suflices to simply 
add a quadratic term to the hamiltonian ( ^.l| ), 



a>l3 



(2.23) 



and for the construction of the reversible mode-coupling 
terms, again all that is required are the following Poisson 
brackets, 

-5"^Ml^'< - S^^M""^ . (2.24) 



Upon inserting ( ^.24 ) into Eq. ( |2.6| ), one readily finds 
the following mode-coupling terms in the equations of 
motion of the order parameter. 



/3 



SH 





(2.25) 



and in the equation of motion for the conserved angular 
momenta, 



SH 
SH 



6H 



SH 

SM^-y 



= -g 



(2.26) 



respectively. Note that as for model J [Eq. ( 2.16 )], here as 
a consequence of the antisymmetry of the Poisson brack- 
ets only the gradient terms in the hamiltonian contribute. 
We remark that in the ordered phase the above reversible 
mode couplings produce propagating Goldstone modes 
with linear dispersion uj{q) oc q. 

Thus we arrive at the following set of coupled nonlinear 
Langevin equations that define the SSS model, 



/5 



dt 



dt 



= -9a 



2o0 



SH[{SS}] , 



with (C"(x,<)) 0, (77"''(x,t)} = 0, and 



(2.27) 

(2.28) 
(2.29) 



(r (x, i)C^(x', t')) = 2Ao<5(x - x')<5(< - t')S 
(x, t)7j^\^', t')) = -2^oV2,5(x - x')<5(i - t') X 



Here we have allowed for violation of detailed balance 
via introducing noise coefficients Aq and Dp which are 
in general not taken equal to AofceTs and D^k^TM, re- 
spectively, where Ts and Tm are the temperatures of the 
heat baths of the order parameter and of the conserved 
variables. Yet we now also have the additional freedom to 
choose the ratio 9o — Ts/Tm different from 1, which cor- 
responds to a violation of the detailed balance in the dy- 
namical coupling of the modes and M^^. We should 
stress again that neither the underlying 0{n) symme- 
try nor the spatial isotropy are affected by this specific 
nonequilibrium perturbation. 

This becomes clear up on co nsidering the Janssen- 
De D omin ic is fu nctional ( 2.14 ) which corresponds to 
Eqs. (2.27)-( p.3C| ). Its harmonic part now reads 



na 



J2 DaMfV^Mf - 

a>f3 

' E (I - DoV^ ) Mo 

a>f3 ^ 



(2.31) 



and can be readily rescaled to the equilibrium form via 
Eq. (2.21) combined with 



Dn 



Dn 



1/2 



a/3 



(2.32) 



Thereby, the relaxation vertex 



a,l3 



attains the new effective coupling 



- Ao 



(2.33) 



(2.34) 



For the mode-coupling terms. 



- fd'^x fdtJ2- 

a,0 



goSoM^'^S^ 



(2.35) 



6 



however, which originally have identical couplings go, the 
effect of this rescaling procedure is to generate two differ- 
ent dynamica l cou pling constants in the first and second 
terms of Eq. (2.35), respectively, namely 



- Do' 



1/2 



90 



9o = ©0 go , 



where 



_>^o Do 

^"-ToD-o 



(2.36) 



(2.37) 



Thus, even if both equations ( 2.27 ) and ( 2.28| ) obey 
detailed balance separately, two different dynamic cou- 
plings will be generated as long as Tg 7^ Tm, and then 
the new variable 80 = Ts/Tm describes the deviation 
from equilibrium. With the two independent couplings 
^0 and g'o, the renormalization group equations will be- 
come different as compared to the equilibrium situation, 
and new critical behavior may be expected at least in 
the extreme cases where the temperature ratio is either 
00 = or 00 = 00. In the following Sec. HI, we shall 
proceed with a detailed investigation of the one-loop flow 
equations of the none qu ilibri um SS S mo del, as given by 
the field theory ( 2.31 ) , (2.33 ), an d (2.35 ), as function of 
the couplings (^.34D, (|2.36|) and (|2.37|). 



III. RENORMALIZATION OF THE 
NONEQUILIBRIUM SSS MODEL 

A. Response functions and Ward identities 

By addi ng sou rce terms to the Janssen-De Dominicis 
functional (2.14), one arrives at the generating functional 



X exp / d^x / dt (h""^^ + /i"Vo ) : (3-1) 



and the (7VA'^)-point correlation functions (cumulants) 
G'^ r A , w can be obtained from In Z via functional deriva- 

fives with respect to the sources ft." and /i" , and then tak- 
ing all h" = h"" = 0. Following the usual field-theoretic 
techniques [ p3|p^ , we furthermore define the generating 
functional for the one-particle irreducible vertex func- 
tions using (pQ — 6 In Z/5h" and ip" = S\nZ/Sh" via the 
Legendre transform 



r[m,{cf>o}]=-^^z[{h'^},m] 



+ 



dt 



(ft"^^-Fft>^) ; (3.2) 



the (A'^A'^)-point vertex functions rp^jv^w then follow via 



functional derivatives of (|3.2D with respect to (f>o and (po . 



With {ip^) = we can write (V'^(x, t)V;^(x', i')) = 
G"^ 7 , (x — x',i — t')S"^ , etc., and upon introduc- 

ing the Fourier transform according to — 
J d'^x J dt'ip°'{:si, t)e^^^'^'^~^^\ one finds the following con- 
nections between the two-point correlation and vertex 
functions. 



G'o^^(q, w) 



(3.3) 
(3.4) 



For the A^iV-point functions with NN > 2, relations 
similar to (3.4) hold, see Eq. (B.9) below. 

In order to assign a meaning to the auxiliary fields, we 
compute the response functions for the SSS mod el by first 
adding external fields to the hamiltonian ( 2.23 ) ||l2| , ^ . 



a a>l3 



, (3.5) 



which produces the following additional terms in the dy- 
namic functional. 



J ^ J- 



d^x / dt 



Xo Y^'^Sg -DoJ2 H'^^'^^K'^ 



a, 13 
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(3.6) 



Therefore the dynamic order parameter susceptibility be- 
comes 



Xo(x-x',t-t')^"^ = ^^^^ 
' 5hP{^',t') 

= Ao(5o"(x,i)^o^(x',0) + 



-goE<'^o(x,t) 



(3.7) 



and similarly the response function for the conserved 
quantities reads 

Xo(x - x', i - t') {5°'-'5P^ - 5°'^5P-<) = 



H-i^=0 



SH'-'^{x',t') 
-Do{AC{x,t)V'M^'{x',t')) 



-2go{Mf{x,t) 



r<7 rtd 



(3.8) 



[Note that X;p(A^o"^(x,i)[M(/''M^'''](x',i')) = 0.] Hence 



we also need cumulants containing composite operators 



^0" = EaiMo'^S^] and Y^"'^ = [<5o"'5'o], as weU as the 
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corresponding vertex functions, which are related to each 
other via 



G'Ssy(q,t^) 



r, 



OSY 



(3.9) 



Using Eqs. (3^) and (3^), we can finaUy write 



(3.10) 



(3.11) 



We conclude this discussion of general properties of the 
SSS model with the derivation of Ward identies which 
are a consequence of the 0{n) symmetry, and the fact 
that the M^^ are the generators of this symmetry group 
| p6| , p^j2l| . As a first version, consider that a spatially 
homogeneous, but time-dependent external field H"^{t) 
is switched on at t = 0. A ccording to Eq. ( |3.5| ) and the 
equation of motion ( 2.27 ), this produces the following 
additional contribution to the expectation value of the 
order parameter component Sq, 

(5o"(x,i))H = -50 fdt'H'^0{t'){Si^{^,t'))^ . (3.12) 



Upon employing suitable variational derivatives, this 
leads to the following relation between the nonlinear sus- 
ceptibility 



Ro S;SA/(x, t; x', t'; x", t") {S^-'d^'' - r^<5^^) 



Shl^{x',t')SH'r^{x",t") 
and the order parameter response 

J d''a;'i?oS;SAf(x,t;O,0;x',i') = 



(3.13) 



-goe{t-t')xo{^,t) . 

(3.14) 



An equivalent Ward identity for vertex functions can 
be obtained by noting that the "mixed" generat- 
ing functional W[{4>^},{<j)^},{H"'f^},{H°''^}], [compare 
Eqs. ( |3.lD , ( ^.2[ )] is invariant with respect to t he fol - 
lowing nontrivial variations corresponding to Eq. (3.12): 

^ ^jjaP^ ^ -eg^Y, H<^0 (t)f^t ||| . HCUCC 



sw 



sw 



, 



(3.15) 



which with jx^^ = S\nZ/ SH"^ translates to a Ward iden- 
tity for the generating functional (3.2) of the vertex func- 
tions. 



2go 



ST 



bo . 



= 0. (3.16) 



fd'^x fdty^ — 



Specifically, this yields 



To M55(q/2, ^/2; q/2, u;/2; -q, -u) = 

= 90 0^)^^OMMi^/^^^/^)^ossi^/'^^^/'^)] ■ (3-17) 



Note that Eqs. (3.14) and (3.17) hold quite independently 



of any detailed-balance condition. 

We can use these Ward identities now to demonstrate 
that the mode-coupling constant g, as a consequence of 
the underlying 0{n) symmetry, does not renormalize [ p^ . 
First, we note that the static response function for the 
conserved angular momenta is exactly 



Xo(q,c^ = 0) = 1 



(3.18) 



as follows from the hamiltonian (2.23) and the fact that 
in the limit w — > there is no coupling between the 
critical fluctuations Sq and Mq^, which is true even 
for our nonequilibrium model. Therefore there cannot 
be any field renormalization for the angular momenta: 
M"l^ = zIPMq'^ (and similarly for M°'f) with Zm = 1. 
Second, as a res ult of the q dependence of the mode- 
coupling vertices ( 2.35 ), to all orders in perturbation the- 
ory 



d{iuj) 



(3.19) 



and hence Z^jZm 



1. We remark that an analogous 
equation for model B leads to the identity z = A — rj 
for the dynamic exponent |12|; a similar result for the 
KPZ equation implies the absence of field renormaliza- 
tions there as we ll |24[|. At last, we utilize the above Ward 
identities (3.14), ([3.17), both of which imply that the 



renormalization factor for the mode-coupling constant is 
identical to Zg = Zm = 1- Physically, this means that 
the reversible mode couplings are not affected by critical 
fluctuations. This fact will lead to certain general iden- 
tities for the dynamic exponent in the scaling regimes, 
see Sec. 



IIIC 



Again, similar Ward identities may be 
derived and corresponding conclusions can be drawn for 
the mode-coupling constant in model J, as mentioned 
above, leading to the exact result z = (d -|- 2 — ry)/2 [ p^ , 
and also for the nonlinearity in the KPZ problem, there 
originating in the Galilean invariance of the equivalent 
Burgers equation, and implying the nontrivial scaling re- 
lation z + X = 2 between the dynamic and roughness 
exponents p4[. 
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B. Renormalization to one— loop order 

Bearing the results of the previous subsection in mind, 
we introduce multiphcatively renormalized fields and pa- 
rameters according to 



~ S ' ■S' ' 

A = Z-^Zj^Ao , D = Zj^Dq , 

A = {ZgZs)-'^''ZxXo , D = ZdDo , 
To = ro - roc , 



T = Z Z.,-Tq^^^ 



u 



Zg^ZuUoAdfi 



d-4 



(3.20) 
(3.21) 
(3.22) 
(3.23) 
(3.24) 



Here, = r(3 - d/2)/2'^-^T:'^^'^ denotes a d-dependent 
geometric factor, and is a momentum scale. Thus 
all the renormalized couplings are dimensionless, as is 

g — go^y^M^'^"'*^^^- Note that both the static and the 
mode-coupling constants uq and go become dimension- 
less at the upper critical dimension dc = 4. We determine 
the renormalization constants (Z factors) by demanding 
that they absorb all the (ultraviolet) divergences in the 
corresponding vertex functions to one-loop order (see the 
Appendix) . We furthermore employ the dimensional reg- 
ularization scheme with minimal subtraction in d = 4 — e 
dimensions, i.e., only include the ultraviolet poles cx 1/e 
in the Z factors, along with their residues in four dimen- 
sions (further details on these procedures can be found 
in Ref. p3| ). In order to avoid the infrared singulari- 
ties near the critical point, we take t — 1 {tq — to 
one-loop order) and q = 0, w = as the normalization 
point. This, of course, follows closely the renormaliza- 
tion procedure for the equilibrium SSS model [T9| ] (see 
also Ref. 

Using rijij(q,a') ~ Z'^^T^ gg{q,u}), the renormaliza- 
tion of the noise strengths Aq and Dq, as well as of the 
diffu sion co nst ant Dq is readily inferred from Eqs. ( [A3| ), 
(A4), and (|A8|), respectively, with the results 



Zd= 1 



e 

2e 
2e 



1 



1 



■Wo 



■wofo , 



Wo fo 6o > 
Wo fo ©0 , 



where we have used the definitions ( 2.37 ) and 



Wo 



\o_ 
Do 



I 2 Do 



(3.25) 
(3.26) 
(3.27) 

(3.28) 



for the ratio of relaxation constants wo and effective dy- 
namical coupling Jo- 

Next, we consider r55(q,w) = {Z ^Zs^'^/'^T^^ gg{c\,u), 
see Eq. ([A5|). First, we determine the fluctuation- 
induced Tc shift roc from the condition of criticality 
Xo(q = 0,LJ = 0)^^ = 0, which because of Eq. (3.1C) 



is equivalent to demanding that gg (Oi 0) 
ro = roc- Eq. (|A^ ) then yields with Eq. (|2.34[) 



for 



■uo 



1 



6 Jk TOc + fc^ 

-(n-l)z(;o/o(l-eo) 



1 



kWo roc + (1 + wo)^^ 



(3.29) 



note th at fo r d < 2 the integrals on the right-hand-side 
of Eq. (3.29) are infrared-divergent, which indicates that 
die = 2 is the lower critical dimension. Evaluating the 
momentum integrals for 2 < d < 4 gives explicitly 



roc 



2Ad 



(d-2)(4-d) 

d/2 



Wo 



1 -I- Wo 



-Uo + (n - 1) X 

2/(4-d) 



fo (1 ' eo) 



(3.30) 



(notice the pole at die = 2 and the essential singularity 
at dc — 4:). The first term here corresponds to the down- 
wards shift of the critical temperature of the <p'^ model; 
the second contribution, which is of purely dynamical ori- 
gin, may either reduce Tc further, namely for Ts < Tm, 
or enhance it with respect to the equilibrium situation, 
iiTs>TM. 

Upon defining tq = tq — toc, th e tru e distance from the 
critical point, and inserting Eq. ( 3.29 ) into (A7), setting 
fQc = + O(uo, go) in the integrals, one finds 



Z\ AitUL 

Zr^ = 1 — 

Zs e 
AdlJ.' 



n - 



A, 



■Uo 



1 



6 e 1 + Wo 

— — -J /o 1 - Bo) • 

(1 + WoY 



Wo fo Oo 
(3.31) 



Then, rendering dV gg{{),Lo) / d{iu}) and 9r^g(q, 0)/9g^ 
finite gives 



iZsZs) 



1/2 



AdfJ- 



[l + wo) 



M/o(i_eo), (3.32) 



Zx = 1 



Ad^J-'' 

e 

AdfJ.' 



n — 1 

1 + Wo 

' in - l)wl 



{i+wor 



Wo fo ©0 - 

fo (1 - ©0) 



(3.33) 



Eq. (3.32) also absorbcs the divergences in the three- 



point function (All), which confirms that indeed Zg 



Eq. (A12) can then be used to determine the still un- 
known field renormalization itself, with the result 



Zs = l- 



Adi^ 
2e 



1 

(1 + wo)' 



■wo/o(l-eo); (3.34) 



notice that Zs ^ ^ and Zg ^ 1 already to one-loop order 
if ©0 7^ 1. At last, the multiplicative renormalization of 
the nonequilibrium SSS model vertex functions is con- 
cluded by rendering the four-point function (A13) finite 
with 
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Zu^= 1 ^ uo H -;— Wo /o 

e 6 e 1 + Wo 

5" Jo i-L - foj - 



e (l + iwo) 



1 + Wo 



uo 



eo(i-eo). (3.35) 



When detailed balance holds, Oq = 1, these one-loop Z 
factors reduce to the well-known equilibrium results [ p^ . 

Whereas the vertex functions and hence also the two- 
point correlation functions are rendered finite with 
the above Z factors, the dynamic response functions 
may require additional additive renormalizations, as a 
consequence of the involved composite operators. The 
r espo nse function for the conserved angular momenta 
(3.11), using Eq. (3.18) for its static limit, can generally 
be written in the following form. 



^o(q,w) 



Ao(q,w)(72 



Ao(q,cj)g2 



(3.36) 



Using Eqs. ([Aq) and (AlO), one finds the one-loop result 



Ao(q,w) = Do 
1 



2 - f fc2 

d 7fcTo + (q/2 + k)2 



1 



■ TO + (q/2 - k)2 -iw/2Ao + To + g74 + 



; (3.37) 



hence, as the ultr aviole t singularity in ( 3.37 ) is absorbed 
by the Z factor (3.27), no additive renormalization is 
needed. This comes as no surprise, as the contribution 
from Eq. ( AlOD is nonsingular. 



However, the integral Eq. (A9) is divergent, and there- 
fore a corresponding additive renormalization has to be 
introduced. The structure of the order parameter sus- 
ceptibility ( |3.10D is 



xo(q,t^) 



Ao(q,w) 



-iw-l- Ao(q,w)/xo(q,0) 



(3.38) 



Here, using Eqs. (A5) and (A9), the static susceptibility 
reads to one-loop order 



Xo(q,0) ^ = To 
n — 1 



n + 2 ^ 
1 TT- uo 



1 



fcfc2(T0 +fc2) 



1 + Wo 



u;^/o(l - Go) X 



k^[woTo - (1 - wo)(q • k) + (1 + wo)((7V4 + fc2)] 

^2 ^ " ~ ^ wo/o(l - eo) X 
1 + Wo 

(1 -wo)(q-k) - (l + wo)gV4 

fcfc2[woTo - (1 - wo)(q-k) + (l + wo)(g2/4 + P)] 

(3.39) 



where Eq. (3.2£) has been inserted, and the renormalized 
Onsager coefficient is 



Ao(q,w) = Ao 



1 + (n - l)wo/o0o 
1 



1 



,To + (q/2 + k)2 

; (3.40) 



-tLo/Do + wo [to + (q/2 + k)2] + (q/2 - k) 



as expecte d, th e abov e multiplicative renormalizations 
with Eqs. (3.31)-(3.34) do not suffice to remove the di- 
vergences in Eqs. ( 3.39| ) and ( 3.40| ). We determine the 
necessary additive renormalization by requiring that 



d 



;^Xo(q, 0)s\L ^ Zs + As , 



and Eqs. ( |3.39| ) and ( |3.34| ) then yield 
AdH"" n - 1 f 1 Wo 



As^ 



(l + Wo)2 



1 + Wo 



(3.41) 



wo/o(l - ©o) 



(3.42) 



Indeed, xo{0,0)-^ and Ao(0,0) are then rendered finite 
with the combinations of Z factors {Zs + As)Zr/Zs and 
Zx/ {Zs + As), respectively. 



C. RG flow equations and fixed points 

The renormalization group equations serve to connect 
the asymptotic theory, where the infrared divergences 
become manifest, with a region in parameter space (in 
our case consisting of {a} = A, D, A, Z?, m, t) where the 
couplings are finite and ordinary "naive" perturbation 
expansion is applicable. They are derived by observing 
that the "bare" vertex functions do not depend on the 
renormalization scale /x, 



d 



'^aS'-M''S-M'(Wo}) - . 



Introducing Wilson's flow functions 



d_ 
9/i 



InZ^ 



In^ 

1 ao 



d_ 



\wZs 



(3.43) 

(3.44) 
(3.45) 



Eq. ( |3.43| ) may be written as a partial differential equa- 
tion for the renormalized vertex functions 



d 



d 



{a} 



(3.46) 
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Note that = Cm = a nd Co = — e /2 as a consequence 
of the exact results in Sec. [II A . Eq. ( 3.46| ) can be solved 
with the method of characteristics /i ii£; this defines 
running couplings as the solutions to the first-order dif- 
ferential RG flow equations 



' dl 



(3.47) 



The solution of the Callan-Symanzik equation ( 3.46| ) 
then reads 



f 1 



exp 



rCs(f) + <5(0 



dl' 



Upon introducing the renormalized ratios 



w — — , B = — — 
AD 

and renormalized effective couplings 
b ^ A 



A 



(3.48) 



(3.49) 



(3.50) 



and collectin g the d efinit ions E qs. ( 3.20)-(3.23) and one- 
loop results (|3^) -(^^ and (|3.31|)-(p.34|)~one finds 



n — 1 



1 



2 (l+w)2 
n — 1 1 — AiW 



(1 + w) 



u;/(l-e) , 

2^/(i-e) , 



Cr = -2 



n + 2 (n-l)z(^ - 

■ u H — 7^ w / (1 — B) 



6 



(3.51) 
(3.52) 
(3.53) 



Cd = - 

Ca = - 



(l + u;)3 
n—\ - n—1 1 — Aw 



1 + w 
2 

71 - 1 



(1 + H' 



(3.55) 



1 + U) 



(l + w)= 



1 



Ci5 = -2^/e 



(3.57) 



Notice that nonzero values of and C^j induce anoma- 
lous noise correlations, while Ca an d determine the 
dynamic critical exponents, see Eq. (3.66) below. 

We furthermore need the flows for the running cou- 
plings v{l), with {w} = w, 0, /, u. 



^^=/3.(^), v(X) = v, 

as given by the beta functions 

d 



= /i 



V ; 



(3.58) 



(3.59) 



with Eqs. ( |3.24D , ( |3.35| ), and (|^)-(|^) these become 
to one-loop order 



/3u.= w (Ca - Cn) 





n 






1- 





e 



(n - l)w^ 



(1-6) 



(3.60) 



/3e= e (C>, - Cz3 - Ca + a 



--w/e(i-e) 



e -f (n - 1) 



1 -I- 7w -I- Aw^ 



(3.61) 



/3j=/(-e + C,)-2CA-CD) 



/ 



2 



1 2 ?i - 1 
- + -^^ 

2 1 + w 



2(n- l)w2 



u'/(l-e) 



(1+U.)3 

n + 8 . 2(n - l)u; 

?y — 

6 (l + w)3 



(3.62) 



-e-f 



«;/(l-e) 



-^^fi-^)W(i-e) 

1 w V u ' ^ ^ 



(3.63) 



We are now ready to explore the fixed points of the 
RG flow e quat io ns, a s given by the zeros of the beta 
functions ( 3.60 )-( 3.63 ). First, we can check that indeed 
for 8* = 1 the equilibrium fixed points (see Ref. [p^) 
emerge. The above flow equations then simplify co nside r- 
ably, and the effective dynamical coupling in Eqs. (3.54)- 
( |3.57 ) becomes 



(3.64) 



because as now = Ca and Cf; = C-D, we can identify 
the noise strenghts A and D with the Onsager coefficients 
A and D, respectively. The corresponding beta function 
for / reads 



/3/ = /(-e-CA-Cz5)=/ 



n- 1 



/ 



(3.65) 



The first equation in (3.65) implies that for any non- 
trivial fixed point < /* < cx) the exact relation 
Ca + Cd = — e = d — 4 holds. Furthermore the analy- 
sis of the RG equation in the vicinity of /* reveals that 
the dynamic exponents for the fluctuations of the order 
parameter and conserved quantities are given by 

zs = 2 + CI , ZM = 2 + Cd , (3.66) 
which then leads to the following identity ]l6| , pT[| 
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Zs + ZM 



(3.67) 



Therefore, in a strong-scaling situation where the char- 
acteristic time scales for the order parameter and angu- 
lar momenta are the same, Q < w* < oo, and hence 
Zs — zm = z^ one finds the well-known exact result 



d/2 



(3.68) 



I ndee d, the above one-loop flow equations ( 3.60| ) and 
(3.65) provide the strong-scaling fixed point 



3 7 /oq — ^ 



(3.69) 



However, in addition there are two nontrivial weak- 
scaling fixed points with 7^ Cdi namely 



with 



and 



Zm— 2 







f* 

J u 



2e 



2n~l 



(3.70) 



2{n - l)e 

2n- 1 
e 



2n- 1 



+ 0(e2) 



00 



/i' = 2e , 



implying that 



Z5 = 2 , ZM = d-2 



(3.71) 



(3.72) 



(3.73) 



Note that at both fixed points the relation zs 4- zm = d 
holds, of course. For the fix ed po int ( 3.72| ) this sum rule 
actually even implies that ( 3.73| ) is probably exact, for 
should vanish if w* = 00. Finally, there is also the 
(Gaussian) model A fixed point with /q = (and w* un- 
specified because now the order parameter and angular 
momenta are d ecoup led), with zs = 2 + 0{e^); yet, ac- 



cording to Eq. ( ^.651 ), it is clearly unstable against / for 
d < 4: £df/di = f3f{e) = -ef at /g* = 0, and hence / will 
increase in the asymptotic limit — > 0. Similarly, both 
weak-scaling fixed points are unstable (to one-loop or- 
der at least) for d < i: Near wj^, = 00 one has (3^ = ew, 
and hence w will decrease as ^ 0, while in the vicin- 
ity of w*, = one finds (3^ = ~£w{2n - 3)/(27i - 1), 
and upon decreasing £, w will go down as well. Sta- 
bility analysis therefore demonstrates that to one-loop 
order the strong-scaling fixed point with z = d/2 is sta- 
ble |jl^,|l8|,^; however, as = is actually close to 
its stability boundary for n = 2, it may well be that to 
higher loop orders the strong-scaling fixed point ac tually 
becomes unstable for the planar model, and in fact (3.6S) 
does not hold |l|,|2|. 

We shall not pursue this issue further here, but rather 
turn to the newly emerging, genuinely nonequilihrium 
fixed points, at which even the static critical behaviour 



might be changed, as opposed to the equilibrium situ- 
ation with O* = 1 for which statics and dynamics de- 
couple, see Eqs. ( 3.51 ), ( |3.53D , and ( 3.63| ). For all the 
above fixed points, we therefore have the nontrivial static 
Heisenberg fixed point (to one-loop order) 



6e 



(3.74) 



which is stable for d < 4, and leads to the critical expo- 
nents of the 0(n)-symmetric c/)"* model 



i/. = -C = 2-(!i±^ 

' ^ I o 



(3.75) 



At the critical point there also appear anomalous long- 
range noise correlations, both for the fluctuations of the 
order parameter and of the conserved fields. For the or- 
der parameter noise, given by r^^(q, lj), one finds 

(r (q, c.)C^(q', Lo')) = A(q, c.)<5(q + q')<5(u; + lu')5'^P , 

(3.76) 

with 

A(q, 0) oc ^ ^(0^ ^) ^ ^(^s-2)/2 ^3 77) 

Similarly, from T £jj;^j{(\,u}) we can infer the noise for the 
generalized angular momenta at the critical point, taking 
into account its diffusive character, 

(7y"'3(q^ c.)r;^*(q', u'))^ A(q, c.)<5(q + c{)5{uj + lo') X 

taking into account its diffusive character, we find the 
limiting behavior 



A(q,0) ocg- 



A(q^0,tJ/g2) oc (tj/g2)^«/2 



(3.79) 



This concludes our discussion of the equilibrium fixed 
points, and we now turn to the two new universal- 
ity classes appearing as a consequence of our genuinely 
nonequilihrium perturbation. 

The b eta function for the temperature ratio Oq = 
Ts/Tm ( 3.61 ) reveals that there can only be nonequilih- 
rium fixed points with either (i) O* = or (ii) 0* = 00, 
which as T5 « Tc effectively either correspond to a 
"renormalized" temperature T^/ = cxd or Tm = 0. In the 
first case, Q* = 0, one finds fi^ = —{n—l)wjw^/ {\+w)^, 
with the s table fixed point w* — 00. Inserting this into 
Eqs. (|3.62D and (|3.63D yields 



e* = 



u*= 2u*H = 



12e 



2(n- 1) ' 



n + 8 



(3.80) 
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The divergence of w* already shows that this fixed point 
describes a kind of weak-scaling be havio r some what sim- 
ilar to the equilibrium fixed poi nt (t3.72 ) wi th (3.73 ). In- 
deed, inserting ( 3.80 ) into Eqs. ( |3.56D and ( 3.57 ) yields 

zs = d/2 , ZM = 2 , (3.81) 



i.e., the dynamic exponent for the order parameter is 
identical with its equilibrium value, while the angular 
momenta are described by mean-field theory (simple dif- 
fusion). Actua lly, a ny no ntrivial fixed point < /* < cxd 
via Eqs. ( |3.62 ) and ( 3.66| ) implies the identity 2zs + zm = 
d -I- 2 -I- Cn , and as fo r Q* = both the anomalous di- 
mensions ( |3.57 ) and ( ^.55 ) should vanish accordi ng to 
the general structure of the couphngs, the result ( 3.81 ) 
is probably exact. There is, however, a nonzero anoma- 
lous dimension for the order parameter noise, ('i = — e; 



therefore the constant Aq in Eq. ( 2.2S| ) is to be replaced 
by a wavevector- and frequency-dependent function, 

(r(q,c.)C^(q',c.')) = A(q,c.)5(q + q')<5(c. + c^')'5"'' , 

(3.82) 

with the singular large-wavelength and low-frequency 
behavior (d < 4) 



A(q,0) ocg''-\ A(0,w) (xw('^-4)/2 , 



(3.83) 



which follows from the matching conditions ii£ ^ q and 
(/i£)^ = cj, respectively. Finally, the new noncquilibrium 
static fixed point u* — 2m^ along with the fact that 
now the d ynamics affects the static anomalous dimen- 
sion (3.53), yield the new static critical exponents 



-G = + O(e2 



1 2(n + 2) 



n + 8 



e + 0(e2). (3.84) 



According to Eq. ( ^.36 ), the fixed point ( 3.8C ) corre- 
sponds to the situation where there exists a coupling 
cx f/o of the order parameter to the angular momenta, 
leading to the dynamic exponent zs = d/2, the genera- 
tion of long-range noise correlations, and even to anoma- 
lous static exponents, yet the dynamics of the conserved 
quantities themselves remains unaffected by the critical 
fluctuations and hence displays mean-field behavior. 

The above analysis has tacitly assumed that w = 
wQ — X/D remains finite for 8* = and w* — oo. 
Howev er, = w{Cx — Co) = —ew, and hence the fixed 
point ( |3.8C| ) is unstable for d < 4 against an increasing 
coupling w, or equivalently, against the generation of the 
new effective coupling 



X^D 



(3.85) 



which characterizes the second noncquilibrium fixed 
point where &* — oo, and which describes a coupling of 



the critical order parameter fluctuations into the diffusive 
dynamics of the angular momenta, but no effect of the 
latter on the equation of motion for the ord er pa ra mete r 
itself. Indeed, the anomalous dimensions ( 3.51 )-( 3.54 ), 
( |3.56 ), ( 3.57 ) all va nish when 8 — > cxd with / held finite, 
and with Eq. (| 3.63p on e finds the standard 4>^ Heisenberg 
static exponents ( |3.7q ) along with model A and purely 
diffusive dynamics for the order parameter and conserved 
quantities, respectively, 



Z5 = 2 + 0(e2), ZM = 2 



(3.86) 



note that w* cannot be specified because of this decou- 
pling of the mo des. The beta function for the effective 
mode couphng ( |3.85 ) becomes f3^ — /(— e -I- 2^^^ — 3Ca — 

Cf)) = /(— e + //2), and hence we arrive at the following 
fixed-point values 



8* 



2e 



''H 



(3.87) 



At this fixed point anomalous noise correlations f or the 
angular momenta emerge, which in analogy with ( 3.82| ) 
can be written in the form 

{r^'^^iq, c.)r;^^(q', u;'))^ A(q, c.)<5(q + q')<5(c^ + iu') X 



repl acing Eg. ( 2.30|). T heir singular behavior follows from 
Eq. (|3!55|) with (|3.87j), 



A(q,0)cx/, p^d-2 



(3.89) 



which is probably an exact result again, because the ex- 
istence of a nontrivial fixed point /* with p — 2 
implies the relation 3zs + p = d + A 
structure of the perturbation theory and the above phys- 



c* 

2q, and both th^ 



ical interpretation require that Cx 



0. Similarly, 



the frequency dependence of A(q, w) displays anomalous 
behavior, but because of the underlying diffusive dynam- 
ics we now have to take the limit q — > more carefully, 
namely with uj /q^ held fixed. One then gets 



A(q^ 0,t^/q2) oc [uj/qy/^ 



(3.90) 



It is a remarkable fact that these anomalous noise correla- 
tions always appear for those degrees of freedom, towards 
which the energy flows, i.e., those quantitites which are 
in contact with the heat bath at lower temperature. It 
shoul d b e not iced that the power laws in Eqs. ( 3.83 ) and 
( |3.89 ), (3.90) are not determined by the corresponding 



dynamic exponents, as opposed to the equilibrium sit- 
uation described by Eqs. (3.77) and ( 3.79D , which are 
consequences of detailed balance. Yet again, inspection 
of Eq. (||6|) in the vicinity of ( ^^ , /3j = -2e/, shows 
that this second new fixed is unstable with respect to the 
flow of /. 

Therefore, as both new noncquilibrium fixed points are 
actually unstable for d < 4, the asymptotic critical behav- 
ior must be governed by the equilibrium strong-scaling 
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fixed point (3.69), (3.74) with dynamic expo nent ( 3.68| ) 
and the usual static critical exponents (|3.75| ). The sta- 
bility against perturbations away from 0* = 1 is readily 
demonstrated by observ ing t hat (3q = — C(l — Q), with 
C > according to Eq. (3.61). Hence any deviation from 
the equilibrium fixed point will be counteracted by the 
flow of the coupling O; i.e., whenever the initial value of O 
is neither zero or infinite, the flow will asymptotically ap- 
proach the stable equilibrium fixed point, and thereby de- 
tailed balance is dynamically restored. This is illustrated 
in Fig. 1^, which displays the location of the fixed points 
and their flow in the space of the dynamical couplings w, 
/, and /; notice that both the equilibrium and "nonequi- 
librium" model A fixed points are represented by lines 
here, because the undetermined value of w* . Concluding 
this section, we remark that at each of the fixed points 
discussed above, including both new nonequilibrium fixed 
points, the anomalous d imen sion stemming from the ad- 
ditive renormalization ( 3.42 ) vanishes, as does Cs from 
the field renormalization. We therefore did not have to 
take its effects into account explicitly. 



IV. SUMMARY AND CONCLUSIONS 



We have studied the critical dynamics of 0{n)- 
symmetric systems, including reversible mode-coupling 
terms, in the framework of effective Langevin equations 
and dynamic field theory, and generalized the equations 
of motions to nonequilibri um s ituations where detailed 
balance is broken. In Sec. |ll B| we have argued that for 
the dynamics of isotropic Heisenberg ferromagnets the 
effect of violating the Einstein relation between the spin 
diffusion constant and the Langevin noise strength can 
be absorbed via rescaling the static nonlinearity and the 
mode-coupling constant. Hence universal properties can- 
not be affected by this specific form of detailed-balance 
violation, and the critical point is described by the usual 
Heisenberg model exponents with the equilibrium dy- 
namic critical exponent z = (d + 2 — ri)/2 jis). This 
fact generalizes previous results regarding the stability 
of the relaxational models A and B against nonequilib- 
rium perturbations |2|j^ to a situation where reversible 
mode-coupling terms are present as well. 

On the other hand, for the critical dynamics of pla- 
nar ferromagnets (model E) [l^ or isotropic antiferro- 
magnets (model G) ||l8|, b oth incorporated in the 0{n)- 
symmetric SSS model ||l5| , such a simple rescaling does 
not remove the effects of detailed-balance violation com- 



pletely, as discussed in Sec. |II Q . This is because there ap- 
pears a new degree of freedom, namely the temperature 
ratio Qq of the heat baths to which the order param- 
eter and the conserved angular momenta are attached. 
This new variable induces different renormalizations for 
the noise strengths as compared to the Onsager coeffi- 
cients (Sec. HIB), and therefore genuinely new dynamic 
and static critical behavior may emerge. Indeed, there 



appear two new fixed points of the resulting RG flow 
equations, describing continuous phase transitions of en- 
tirely nonequilibrium character, namely corresponding to 
either G = or O = oo, where in the former case even 
the static critical exponents become modified. The en- 
suing dynamic exponents may be interpreted physically 
by noting that in both cases either the coupling of the 
order parameter into the diffusion equation for the con- 
served fields vanishes, while the reverse coupling remains 
effective, or vice versa. However, the stability analy- 



sis in Sec. HI C shows that both of these fixed points 



are actually unstable, and provided < Go < oo the 
asymptotic critical behavior is governed by the nontrivial 
strong-scaling equilibrium fixed point of the SSS model 
characterized again by the static exponents of the 0{n)- 
symmetric Heisenberg model and by the dynamic expo- 
nent z = d/2 ||,|l9|. 

This suggests that the role of detailed balance is actu- 
ally a weaker one as compared to internal symmetries; for 
while breaking a discrete or continuous symmetry typi- 
cally results in a change of the universality class, detailed 
balance becomes restored here at the critical point, ob- 
viously as a consequence of the underlying 0(n) rotation 
symmetry and the spatial isotropy of the model. We 
remark that there are some notable exceptions, though, 
for the relevance of symmetry-breaking terms. E.g., the 
0{n) symmetry of the Heisenberg model ( |2.lD is restored 
as a consequence of the large critical fluctuations, even 
when cubic anisotropics are added, if n < ric (with rtc = 4 
to one-loop order) |2^. We should also cautiously state 
that our above results relied on the one-loop approxima- 
tion only, and specifically the stability boundaries might 
change in more accurate calculations, as seems indeed to 
be the case for the equilibrium planar model p2[| . How- 
ever, in the beta functions Pw and /3 f describing the in- 
stabilities of the fixed points ( 3.80| ) and (3.87), respec- 
tively, no dangerous n-dependences appeared, and it is 
therefore probably safe to say that these nonequilibrium 
fixed points will remain unstable even to higher orders 
in perturbation theory, albeit for n = 2 the asymptotic 
region may ultimately be governed by the weak-scaling 
equilibrium fixed point ( |3.70 ), see also Fig. |l|. 

The remarkable result that violation of detailed bal- 
ance appears to be an irrelevant perturbation in the RG 
sense, at least in rotation-invariant and isotropic sys- 
tems, is of course strengthening the notion of universal- 
ity even in nonequilibrium situations. Probably in many 
experiments probing the critical dynamics with electro- 
magnetic radiation or neutron scattering, the system is 
not perfectly thermalized and variations in the effective 
temperatures for the different degrees of freedom cannot 
be avoided completely due to the inevitable critical slow- 
ing down which prevents fast relaxation processes. Now 
if such a perturbation were relevant, its effects would 
become enhanced drastically in the vicinity of the phase 
transition, however small the initial deviations from equi- 
librium might have been. Such a behavior is apparently 
not observed; however, it might be interesting to prepare 
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a non-thermalized system on purpose, say by arranging 
for the in-plane spins and in a planar ferromagnet 
to be on a different temperature than S^, conceivably 
attainable with polarized electromagnetic radiation. An- 
other possibility would be to introduce a long wavelength 
magnetic field with random time variation. In an anti- 
ferromagnetic material, its effect would cancel for the 
staggered magnetization but the coupling to the magne- 
tization would mimic a situation with Tm — oo. For both 
experimental realizations, or at least in equivalent com- 
puter simulations studies, the effects of the above two 
new fixed points may then be explored, for these should 
manifest themselves in nontrivial crossover behavior, and 
perhaps even the associated anomalous critical exponents 
and noise correlations might be detected. 

We finally remark again that the specific violation of 
detailed balance investigated here was isotropic in char- 
acter, thereby disturbing neither of the underlying sym- 
metries of the SSS model, namely 0{n) symmetry and 
its spatial isotropy. Obviously, it would be interesting to 
see if the above stability against nonequilibrium pertur- 
bations persists even when the detailed-balance violation 



is applied in an anisotropic manner, e.g., by coupling the 
order parameter to conserved angular momenta as above, 
but arranging the noise strength of the conserved fields 
to be related to different temperatures in different space 
directions. The study of this, or similar direction- or 
scale-dependent nonequilibrium perturbations provides 
a promising venue for further research. 
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APPENDIX A: EXPLICIT ONE-LOOP RESULTS FOR THE VERTEX FUNCTIONS 



In this Appendix, we present a list of the results to one-loop order in perturbation theory for those vertex func- 
tion which are required for the r enorm a lizati on of t he no nequilibrium SSS model defined by Eqs. ( 2.27[ )- (|2.3(]| ), or, 
equivalently, by the field theory ( ^.31 ), (2.33), and ( ^.35 ). In all the subsequent expressions, the internal frequency 
integrations have already been carried out via the residue theorem. We use the abbreviation J^, . . . = (27r)^'' / d'^k .... 
We do not explicitly provide the Feynman diagrams themselves, as they are identical with those of the equilibrium 
model |l|]. 

We begin with the two-point vertex functions renormalizing the noise strengths Aq and Dq, 



1 + (n - l)gl 



Dn 



XoDo 7fero + (q/2-fk)2 
i?o(q/2-k)2 



Ao[ro + (q/2 + k)2] 



{lu - iXo[ro + (q/2 + k)2])2 + Dl{q/2 - k)4 [to + zi?o(q/2 - k)2]2 + A2[ro + (q/2 + k) 



- MM 



2 ^0 



Ao-Dq Jk 

1 



(q • k)^ X 



1 



212 



ro + (q/2 + k)2 {uj - i\o[ro + (q/2 + k)2])2 + A2[ro + (q/2 - k) 

1 1 
+ ro + (q/2 - k)2 {oj + zAo[ro + (q/2 - k)2])2 + A2[ro + (q/2 + k)2]2 



Consequently, we have 

ros5(0>0) = -2Ao 
d 



gq2 a MM 



(q,o) 



q=0 



1 + in- 1)50^ 
= -2Do 1 



Dn 



1 



1 



Ao-Do Jk ro + Ao(ro + k^) + Dok^ 
2 2 Aq f k^ 



-,9o- 



XlDo Jk (ro + k^f 



(Al) 



(A2) 



(A3) 
(A4) 
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For the computation of the response functions, one needs the two-point functions 



6 



+ (n l)(ro + q^)gl ^ ^^^^ ^ _^ ^^^^^ _^ ^^^^ + k)2] + Doiq/2 - k)2 

, / 1 \ 2 ^0 / , ^0 Do 



1 



1 



Do \ XoDo Jk + AoN + (q/2 + k)^] + Do{q,/2 - k)2 ' 



roMM(q>-) -^^ + Doq' - ^ (q/2 + k)2 + 2Ao(ro +qyA + k^)' 



(q.k) 



1 



(A5) 
(A6) 



specifically, 

ro55(O'0) = Ao 



ro(l + in-l)g'o^ 



Jkra + Ao(ro + fc2) + DoA:^ 



71 + 2 Ao 

7— Mo 



6 Ao Jk ro + k"^ 



Aoi'o V ^obo Jk + fc^) + Dok^ 



d_ 
dq 



a„2 ^0 Mm('I'O) 



q=0 



Ao 



d"" XlDo Jk (ro + k^r 



Furthermore, the following vertex functions containing composite operators are required [note Eq. (|3.9| 



"os[MS](q'^) = -('^- 1)30 
Ao 



roM[ss](q'^) - 2.90- 



Ao f 1 1 

Ao Jkro + (q/2 + k)2 -ico + Xo[ro + (q/2 + k)2] + Do{q/2 - k)2 ' 
(q.k) 1 



Ao Jk ro + (q/2 + k)2 -iuj + 2Ao(r-o + q^/4 + fc2) 



(A7) 
(A8) 

(A9) 
(AlO) 



Calculating the three- and four-point functions is already a rather tedious task even to one-loop order, and we 
merely quote the final results needed for renormalization purposes: 



rossM(0:0;0,0;0,0) = -go 



1 - (n - 1)5, 



2D0 f 1 

"l?oyfc[Ao(ro+F) + Z?oA;2]2 



+ 



Jk 



1 



d 

5(q • p) 



OMSS 



(-q,0;q/2-p,0;q/2 + p,0) 



q=p=0 



2.90 



fcro + fc2 [Ao(ro + fc2) + £)ofc2]2 
2, obo 



(All) 



^ d^'' Ao jk ro + k^ [Ao(ro + k^) + Dok^? 



+ 3('^- l)5o 



1 



VDo 

XI Jkiro + k^r[Xoiro + k^) + Dok^Y 



(A12) 



and finally 



ro5sss(0'0;0'0;0'0;0'0) = Aouo 



, n + 8 Ao 
1 — —uo 



6 Ao Vfc(ro + fc2)2 

+ (n-l)9o^fl-^)^/^ 
V Ao-DoMo/ AoDo JfcJ'o + 



+ 



+ fc2 Xo{ro + k^) + Dok^ 



-("-l)5o 



^ I 3.9§ \ £0 r 1 

XqDoUo) Do Jk [Xoiro + k"^] 



+ (n 



fc2) +Dofc2]2 

1 



+ 



ro + fc2 [Ao(ro + fc2)+DoA:2] = 



16 



H Do 



XoDoUo XI Jk {ro + k^)^ Xo{ro + k^) + D^k^ 



-(n-1) 

-(n-1) 



Ik'' 



XoDoUo Xo Aro + k^ [Xo{ro + k^) + Dok'^f 

H ~XoDo f 1 



+ 



XoDoUo XI yfc(ro + A;2)2 [Ao(ro + fc2) + Dofc2]2 



(A13) 
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FIGURE CAPTION: 

FIG. 1. One-loop flow diagram for the nonequilibrium SSS model in the space of dynamical couplings w/{\-\-w), //(I -I- /), 
and //2 (displayed for the case n = 2, e = 1). Asymptotically, the equilibrium strong-scaling fixed point w*q = 2n — 3, /*q = e, 
/eq = £/(2n — 3) is stable; in the figure, it is located at the point (1/2, 1/2, 1/2) (full circle). 
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